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Abstract 

The vast majority of solids in Nature are strongly disordered or amorphous. Many of these 
materials show "universality" in their low-temperature elastic properties, despite enormous differ- 
ences in their microscopic structure^. Discovered 40 years ago^, this universality has emerged as 
one of the outstanding unsolved mysteries in condensed matter physics^. Here we show that the 
mere presence of strong disorder and tunneling states symmetric under spatial inversion necessarily 
gives universality. Such tunneling states interact weakly with phonons, yet gap the non-symmetric 
tunneling states above a crossover energy. Thus neither glassiness nor any amorphous structure 
are required for universality, only strong disorder - a result concordant with recent experiments 2 -. 
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Some features of the universality are well known. Thus, for a large variety of disordered 
and amorphous solids, for temperature T < Tjj ~ 3K, the specific heat C V (T) oc T a with 
a ~ 1, the thermal conductivity k(T) oc with (3 ~ 2, and the internal friction or 
"quality factor" Q(T) = 2irl(T)/ A(T), dividing the phonon mean free path l(T) by the 
phonon thermal wavelength A(T), goes to a constant "plateau" value Q Q ~ 10 3 . Both 
T[/ and Q vary somewhat between materials (by up to 2 — 3 on either side of the quoted 
values). However, the similarity of a, /3, Tu, and Q Q in materials whose microscopic structure 
ranges between impurities in ordered crystals to completely disordered amorphous glasses, 
is really striking, and suggestive of a fundamental mechanism dictating the low temperature 
characteristics of disordered solids. Note that above Tjj things change rapidly - one finds a 
hump in C V (T), a plateau in k(T), and Q(T) falls very dramatically to Q/2n = //A ~ 0(1) 
(see, eg., ref.-). 

In this Letter we present a novel approach to the problem, relating the universal properties 
to the symmetry of tunneling states under local inversion. Tunneling states symmetric under 
inversion have intrinsically weak interaction with the phonon field, defining the relevant 
universal small parameter as the ratio between the elastic and Coulomb energies in the solid. 
Universality then arises from the generality of the symmetry argument, the robustness of the 
energy scales involved, and a generic gapping of the asymmetric tunneling states resulting 
from inter-defect correlations. Being a local mechanism, our theory may shed light also on 
the question of the microscopic structure of amorphous solids. 

Theoretical analyses of the universality began with the influential " Standard Tunneling" 
("ST") model^, in which the low-energy modes are described by a set of non-interacting lo- 
calised two-level systems (TLSs), which however interact weakly with phonons. One assumes 
a flat broad distribution of bias energies, with a probability distribution P(e, A G ) ~ P /A. 
over tunneling amplitude A Q and bias energy e. The central dimensionless parameter of 
the ST model is the "tunneling strength" C Q = P {^ 2 '/ 'pc 2 ), where 7 is the defect-phonon 
coupling, p the mass density, and c the phonon velocity, this model gives C V (T) = a Q T, 
k(T) = (3 T 2 with (3 oc 1/C , and Q = 2/nC , so that A/7 = ir 2 C ; and l(T), A(T) are both 
oc 1/T. Thus C Q dictates phonon attenuation, internal friction, and thermal conductivity. 

Experimentally, C Q ~ 10~ 3 in most materials, even though pc 2 varies by four orders of 
magnitude. There is no reason for this remarkable result to hold in the ST model, and 
for this reason most later theories^^ 1 ^^ 1 ^^ 3 -, recognising the importance of inter-defect 
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interactions, have tried to explain universality as a consequence of these. However this 
leads to a further problem, because in such theories the density of states n(E) = pc 2 /'-/ 2 = 
P /C . Typically P Q ~ O.ln(E)^^: This then implies C Q ~ 0.1, a huge discrepancy with 
experiment. Such problems have led some to argue^i&ii that the defect modes may not be 
TLSs at all, but something quite different. 

• : • 
••• 

(a) (b) 

FIG. 1: Structure of crystals doped with impurities, (a) The 3-dimensional structure of the alkali 
halide KBr with two CN impurities, which align along one of the 3 crystal axes (distorting the 
nearby lattice), (b) A 2-dimensional system in which a small impurity can occupy one of the 
4 positions shown by red and yellow circles; states of the same colour are related by inversion 
symmetry. The low-energy states of this 4-state system are isomorphic to those of an orientational 
impurity like CN which can rotate 90° between a "red orientation" and a yellow one, and flip 
through 180°, whilst keeping the same colour. States of the same colour form a "r-pair", and 
transitions between the two states of the j'-th r-pair are generated by the operator fj. Transitions 
between states having two different colours are mediated by the 5-spin operators. 

Almost all these approaches have assumed glassy behaviour and/or interactions must 
be associated with universality. However amorphous systems without universality have 
been founcU^, and experiments in ion-implanted crystalline Si showed that universality was 
unrelated to the degree of amorphicity^. Universal properties have also been observed in 
the ferroelectric phase^ 1 ^, and recent evidence suggests that strong random strain fields 
acting on tunneling defects, rather than interactions between the defects, are responsible for 
universality^^^i. This is most obvious for impurities in otherwise crystalline hosts. Thus, 
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FIG. 2: Interaction of S 1 - and r-spins with elastic strain. In (a) we show the lattice distorted by a 
passing phonon; this breaks the degeneracy between states of different colours (ie., it couples to the 
S-spins), but not between states of the same color (ie., it does not couple directly to the r-spins, 
to first order in the gradient of the lattice displacement). In (b) and (c) we see how the degeneracy 
between 2 states of a r-pair is broken either in a system in which there is no long-range lattice 
structure (ie., an amorphous system) or by the presence of a nearby impurity in a lattice. In this 
case the r-pair does interact, albeit weakly, with a strain field (shown here as a passing phonon). 

in pure KBr doped with CN impurities, universality requires very high CN concentration 
x (x > 0.2); the impurities then create sufficiently strong random strain fields. However CN 
impurities in (KBr) y (KCl)i- y (where strong random strains already exist) give universality 
for very small x; but now Q Q is x dependent, and for i < 1 one finds Q Q ^> 10 3 . Thus one 
can argue^ that experiments on universality should focus on disordered crystals or crystals 
with impurities, rather than on amorphous systems. 
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FIG. 3: Densities of states of r and S TLSs. General behaviour of the mean field densities of 
states rig(E), n°(E), and the true densities of states ns(E), n T (E), once inter-spin correlations 
are included. These results are shown as functions of energy (measured in temperature units) for 
a typical value of J D ~ 300 K; both axes are plotted logarithmically. The sharp gap-like feature 
in ns{E) appears at a temperature Tjj ~ g,J , where the r-spin levels first appear. Tjj is to be 
identified with the crossover temperature to universality. n T (E) ~ n°{E) is nearly unchanged by 
inter-spin correlations, see text. 

We now outline our theory for disordered crystals, and without loss of generality we 
consider the well studiecU^^yLi^ eX ample of KBr : CN. The applicability of our theory 
to amorphous solids is discussed at the end of the Letter. Fig. 1 shows KBr : CN, with 
the CN impurities oriented along the 3 KBr crystal axes. Each CN then has 6 degenerate 
states, which divide into 3 "r pairs" , of states related to each other by inversion. Transitions 
(90° rotations) between different pairs are mediated by a pseudospin S. For simplicity we 
consider r and S to be spin-half Pauli operators (see Supplementary Information). 

S*-spins differ crucially from r-spins in their interaction with an elastic strain field u Q/ g(r). 
An isolated Sj interacts directly with this field, whereas in a pure lattice an isolated Tj 
can interact only with the gradient of the strain^: inversion symmetry rules out any direct 
interaction (see Fig. 2). However disorder breaks the exact inversion symmetry of the r- 
spins, inducing a weak direct interaction with the strain field. The TLS-phonon interaction 
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(a) (b) 

FIG. 4: The true S'-spin density of states ns(E), as a function of energy for different short-distance 
behaviours of the interaction. We normalise ns(E), defining P(E) = ns(E)/ns(E), where E is an 
energy a few times larger than Tjj; essentially ns(E) is the S'-spin density of states obtained if one 
includes S — S correlations but ignores S — r correlations, so P(E) measures the effect of S — r 
correlations. We use the interaction form = Ci^gJo/iRfk + &o) discussed in the text, where Cj/% 
is a random variable and a Q is a short distance cutoff parameter. Then curves for P(E) are shown 
for different a Q (in units of a), where E is measured in units of E T , the calculated average energy 
for r excitations. Tjj is defined by the energy at which n T (E)^ = ns(E)^, i.e. by P{E) ~ 5g. 
Impurity concentration is x = 0.5 and lattice size is 13 3 cells (giving ~ 4400 TLSs). Similar results 
are obtained for x = 0.2 (see Supplementary Information). For both concentrations we find that 
T\j « 0.2E T for 1 < a Q < 6, i.e. Ty ~ gj . 

then takes the form (see Supplementary Information) 

V» = EE W? S i + + vS aP ] u aP {v 3 ) (1) 

3 a, /3 

Here 7 S w Ec, where Eq ~ 5eV is the typical Coulomb energy in a solid; and the "volume 
energy distortion" parameter t] is comparable to, and usually somewhat smaller than 7 S . 
The parameter 7 w ( r j) varies from site to site, and depends on strain disorder. We define 
a parameter 7 W , the mean variance of 7 w (i\,); it increases with disorder from zero for a 
crystalline host, saturating for a strongly-disordered strain field to 7w ~ E$, where E$ ~ 
O.leV is the typical elastic energy in a solid. Our estimates for 7 S and 7 W are further 
discussed in the Supplementary Information, and are in agreement with recent ab-initio and 
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DFT calculations on KBr:CN^. The small dimensionless parameter g = 7 w /7s = E^/Ec ~ 
(1 — 3) x 10 -2 measures the deviation from local inversion symmetry caused by disorder, 
and plays a crucial role in what follows. Note that r-spins scatter phonons very weakly 
compared to S-spins, down by a factor g 2 ~ 10~ 4 — 10 -3 . 

The defect-phonon interactions immediately generate phonon-mediated interactions be- 
tween all the spins^ 4 - 1 ^, along with random static strain fields acting on them^i. The S — t 
system then has an effective interaction Hamiltonian of form 

H sr = }2-'u Ss f s i + J ^ s * r ! + ( 2 ) 

plus a random field term Hrf = ^2j[hj +fijT*]. We now show how the above Hamiltonian 
necessarily gives the smallness and universality of the tunneling strength C Q . 

The interaction terms J!jf, with a, b = 5 or r, can be written in the form ~ 
Jo b Ctjf(Rij/a), where the ~ 0(1) are complicated functions of angle^ 4 - 1 ^, and the radial 
functions f(Rij) depend on distance R^ like (a/Rij) 3 when Rij 3> a, (a being the lattice 
constant). The three characteristic interaction energy scales J^ S ,J^ T , and J T T are related 
to each other by J T T ~ gJ^ T ~ g 2 J^ s , implying a single fundamental interaction scale 
J = Jq S . Now J = 7 2 /pc 2 i?Q ~ 300 K, (where Rq is the typical distance between nearby 
TLSs), and it does not vary dramatically between different solids. This experimental fact 2 - 1 ^ 4 - 
is also part of the universality picture. The random field strengths are also governed by J Q 
and g\ their typical strengths are hj ~ J and ~ gJ<r^- 

In a naive mean-field analysis the effect of the interactions in equation ([2]) and the ran- 
dom fields is to spread the S'-spin levels and r-spin levels over energy ranges J a and gJ D , 
respectively. One then gets mean field densities of states n° s (E), n°(E) of Gaussian form, 
with rig ~ 1/J Rq for E < J D , and n° ~ l/gJ Rl for E < gJ Q (see Fig. 3). Note that even 
though n°(E) ^> n° s (E) for E < gJ , by a factor 1/g, the S-spins still dominate the phonon 
scattering, because their scattering rate is far higher (by a factor 1/g 2 )- Thus in mean field 
theory, the S-spins dominate the acoustic properties at all energies. 

However this mean field picture is radically altered by inter-defect correlations. We 
summarise the changes as follows (see Fig. 3): (i) the true S'-spin density of states ris(E) 
now shows a slow fall-off from the mean-field result for gJ Q < E < J Q ; then, for E < gJ Q ~ 
3 — 10/^, it falls precipitously, so that for energies <C gJ a , the S'-spin states have essentially 
disappeared, (ii) The true r-spin density of states n T (E) is hardly altered - it does show 
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a weak dip, but only below the much lower energy J T T ~ g 2 J a ~ O.IK, (iii) Finally, the 
phonon spectrum is only weakly altered. 

The energy gj now acquires a physical meaning - it is the energy below which the r-spin 
levels appear, and the S'-spin levels disappear. This rapid crossover is responsible for the 
abrupt change in the system properties at a temperature ~ gJ Q , which we will now identify 
with Ty. The mechanism governing the catastrophic reduction in n,s(E) below Ty is a more 
elaborate version of the "Coulomb gap" mechanism which exists in disordered conductors^ 
(adapted here to the case of 2 different spins and to dipolar interactions). Level repulsion 
caused by the inter-spin correlations forces states away from low energy - the strength of 
this effect depends both on the strength of the interactions, and the density of states of 
the levels. The effect on n T (E) is extremely small, because there are very few \ow-E S- 
states, and because the r — r interactions are so weak (see Supplementary Information for a 
rigorous derivation). The S — S interactions are stronger and are responsible for the slow fall 
in ris(E), but the r-spin levels have a huge effect on the rather small number of low-energy 
S'-levels because their density is so high. 

To evaluate ns(E) we assume a form J^ r = gcik-Jo/ {.R\ k + a 3 ,), which accounts for the 
elastic interactions at short distances (other short-distance forms produce similar results). 
We define Ty by the condition ^\n r {Ty) = j 2 ns(Ty), dictating the crossover between r 
dominated and 5* dominated phonon attenuation. The key result of this calculation (see 
Fig. 4 and Supplementary Information) is that although the shape of n(S) depends on 
a Q , the universal regime exists for any value of the length scale a Q , and Ty is only weakly 
dependent on it for typical values a Q /a ~ 1 — 6; it is ~ 0.2(E r ) ~ gJ Q , where (E T ) is the 
calculated average excitation energy of the r TLSs. 

Thus, below Ty one only sees a set of r-spins, weakly coupled to each other and to 
phonons. Above Ty the S'-spins dominate, along with a plentiful supply of phonons, to 
which they couple strongly. Only in the crossover regime around Ty do S — r interactions 
play any role; and although their detailed form influences the shape of this crossover, it does 
not affect the physics above and below Ty. 

Below Ty the system of r-spins behaves essentially as a set of non-interacting tunneling 
states, which by symmetry are intrinsically TLSs. Thus we may now analyse this system 
using the parameters defined in the ST model. Having n T (E) = pc 2 /^^ and P a ~ O.ln(E) 
we obtain C T Q ~ 0.1n T (E)^/ pc 2 ~ 0.l7 w /7 s ~ O.lg ~ 10 -3 . In comparison to S (asymmet- 
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ric, or standard) TLSs, phonon scattering by each r TLS is smaller by a factor of g 2 , and 
their DOS is larger by a factor of l/g, resulting in C T a ~ gC^. Thus we recover universality, 
now with the correct experimental value for C a . 

The universality then occurs because neither the parameter g = Eq,/Ec nor the energy 
scale J vary much between different strongly-disordered materials. Both the crossover 
energy Ty ~ gJ Q and the parameter C T determining the universal properties must then be 
very similar for most systems; their small size comes from the small value of g. We stress 
that glassiness, and the mechanisms leading to it, are clearly not required in our theory. 
If glassiness does occur, it has no relation to universality here except insofar as the glass 
temperature T G ~ J Q (so that Tjj ~ gTo). 

Finally, let us go beyond the model of a lattice system distorted by impurities, including 
some features omitted above. 

(i) Mixed Crystals: Here the r-impurity concentration x <C x c , but the host is now 
strongly disordered (eg., in KBr:KCl with CN impurities). Thus strain disorder is strong even 
when x = 0; it comes now from the random field exerted by the random host atoms, operating 
through the volume term in equation (CTJ, rather than from interactions with S impurities. 
The model then predicts a reduction in n T by a factor x/x c , so that C Q ~ 0.1gx/x c , in 
agreement with experiment. Conversely, in dilute single crystals 7 W is typically small (oc x 4//3 , 
see Supplementary Information), has a broad distribution, and is correlated with the TLS 
energy: then universality is lost. 

(ii) Electric Dipole Interactions: Both r and 5* defects will have electric dipole moments, 
leading to electric dipole mediated interactions whose strength at nearest impurity distance 
is J ee . Thus J Q > J ee almost always, but the ratio J ee / Jst = Jee/gJo can be small or large, 
and typically J ee > J TT = g 2 J - If Jee < gJo, then all the arguments above go through, but 
now the r — r interaction strength is J ee , and this will change the temperature at which 
the very weak dipole 'depression' appears in n T (E) to J ee . Such dipole depressions in the 
density of states are seen in some experiments^ below 1 K. On the other hand if J ee > gJ Q , 
electric dipole interactions also change the S — r interactions. This does not change the 
basic picture, but now one finds T v ~ J ee , and C Q is reduced to C Q ~ g 2 J /Jee- 

(iii) Amorphous systems and glasses: While amorphous materials show no long range 
order, local order with a lattice bond coordination number does exist except at defect sites, 
and inter-site distances are hardly altered from lattice values. This is true even in strongly 
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amorphous Si, where 8a/ a ~ 0.03 for nearest neighbours, and ~ 0.1 for 2nd nearest neigh- 
bours (a is the regular lattice value as before) and the nearest-neighbour bond angle differs 
by a maximum of 10%. Now our results depend essentially on the distinction between S and 
r defects, and rely on local properties of the system (e.g. contributions to 7 W from disorder 
diminish like ~ 1/r 4 ; see Supplementary Information). Thus neither J Q nor g are strongly 
affected by amorphicity. 

This allows us to make some clear predictions for amorphous systems. First, if inversion- 
symmetric TLSs do exist therein, then our model predicts amorphous systems will also show 
universality with similar values of Tjj and C Q . Second, one should find two types of defect 
excitations, differing by ~ 10 2 in their coupling to phonons, and with the densities of states 
given in Figs. 3, 4. Conversely, universality should be absent in amorphous systems without 
inversion-symmetric TLSs. 
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I. SUPPLEMENTARY INFORMATION 

A. Mapping to the system of S and r-spins 

We consider first the 2-dimensional example given in the text, where the orientable CN 
impurity can lie in one of 4 directions in the host lattice. A similar problem is posed by a 
simple point impurity which may take up one of 4 off-centre positions in the host. One may 
now choose 2 pairs of states from these 4 states; members of each pair are related by inversion 
symmetry, and we refer to these as r-spins; they are described by the Standard Tunneling 
(ST) model. For a single site in the 2-dimensional host, there are two r-spins, r„, with Hilbert 
spaces T u , v — 1,2. Transitions between the 2 states of a given r spin (corresponding to 
180° flips for the CN impurity, or tunneling between 2 positions on opposite sides of the 
centre of symmetry, for simple impurities) are mediated by f x . There are also transitions 
which take the system between T x and T 2 (these correspond to 90° rotations for the CN 
impurities, and hopping between nearest-neighbour off-centre points for a simple impurity). 
Now if we ask what is the interaction between the defect and a phonon field, we find that it 
must have the form 



12 



rf (r>f + Catrtefau? r* (SI) 

where n = 1 . . .1 enumerates the pairs of states (for a 3-d cubic alkali halide lattice, there are 
I = 3 r spins per site). As noted in the text, the strain field couples directly to transitions 
between the subspaces, but only its gradient can couple to individual r variables, because 
of the inversion symmetry. For n = 2 we can write 

rf (r,-) - T a J{v 3 ) = -yf[S}-%]. (S2) 

The resulting analogy to a system of 2 spin half particles is not exact. For example, the 
operator S'Jt? would flip the spin through 180° and then rotate it through 90° in one direc- 
tion, but this is then equivalent to the rotation by 90° in the opposite direction. However, 
in this paper we are only interested in evaluating the effective interactions between defect 
degrees of freedom, for which we only need the Ising variables Sj, f?. Also, the reduction of 
the larger n space of the S degree of freedom to a spin half TLS does not affect the physics 
below T[/, and introduces only small quantitative change to Tjj itself, because of the sharp 
energy dependence of N$(E) at the relevant low energies (see Fig. 4 in the main text). 



B. TLS-phonon interactions and the random strain field 

We wish to show that in the presence of disorder, r-TLSs acquire a finite effective inter- 
action with the first derivative of the phonon displacement, given by 

where u a p{rj) = {dXj a /dxjp) is the strain tensor, defined in the usual way as the gradient at 
site j of the lattice position vector Xj. To derive this we consider the interaction Hamiltonian 
between phonons and defects, of which the relevant terms are 



H in t = 



(S4) 

The first term is a volume distortion term. The second term denotes the interaction of the 
S degree of freedom with the phonon strain, and is just that discussed above (see Eq. (1521 ); 
both terms have been discussed exhaustively in the literature for crystals (see e.g. Refs.- 1 ^). 
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FIG. SI: Interaction of the symmetric TLSs with the gradient of the strain. The degeneracy of 
states related by inversion symmetry is broken by the second derivative of the phonon displacement. 

The third term denotes the interaction of the r-spins with the gradient of the phonon strain 
(see Fig. SI); it was derived for a crystal in ref.-. The last term, 2nd-order in the strain, 
includes the modification of the elastic modulus by the impurities, and is given in its general 
form. In this microscopic Hamiltonian there is no direct interaction between the r-spins and 
u® 13 ; as noted in the text and in ref.-, this is forbidden by symmetry in a pure crystal. 

A simple way to estimate the effect of strain disorder on the interaction between S and r 
impurities, is to re-derive this interaction with a 3rd impurity introduced into the calculation. 
Without this 3rd impurity, the effective TfS? between two distant impurities is ~ 1/-R 4 , in 
2nd-order perturbation theory^. We now calculate the effective t^S* interaction in third 
order in H mt , and deduce 7 W from the coefficient of the term behaving as 1/R 3 for large 
distances^. We assume the 3rd impurity k is much closer to impurity i, to model deviation 
from inversion symmetry resulting from strong disorder strain. 

We write the Fourier transform of the lattice displacement field as 



where e ?)/iiQ , is a phonon polarization index. We then minimise the total potential energy, 
i.e. the sum of the bare phonon potential energy plus the interaction term above, to find 




(S5) 
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the resulting distortion in the lattice^. This distortion is then found to be 

SX e (x) = \ ^(cDf^e^e^ + H.c), (S6) 



2 

qfj, 



where we have 



SX W = - ( £ C^e^qsq^r* + i ife^qpe-^S^A . (S7) 

7,(5,r; a/3 / 



The effect of this distortion on the 3rd impurity k comes only from the fourth term in 
H int , which can be written as 



dx k(j , dx k ^ 

P<PXW 



where 
d5X kp 



dx 



NMu 2 ( ^ ^i SVe i^ s ^ sin <l( x k ~ Xi)7f + Y if^aQp cos q(x k - Xj)S? J B w q^. 



(S9) 

We now evaluate the term in Eq. flS8l) proportional to rf5*|. Defining 

^ = £C^VrW-? (S10) 

-ySr) 

D w = lf e w<* <lP s h ( s 1 1 ) 

a/3 

we find the change in the t^S*- interaction caused by the impurity at site k to be 

E k (i,j) = 2J2 C PM( r k) NMtf Aqp Sin ^ Xk ~ S Dq,f1 ' C ° S ~ 

(S12) 

In the acoustic approximation one finds^ 

= ^ < S13 » 

where ~ 1 is a complicated function of angle, which depends on the position of the 
impurities and takes either sign. It can then be treated as a random variable. Here the 
acoustic approximation is not valid, because Ri k ~ a. However, we are only interested here 
in an order of magnitude, and using the fact that Rn- <C Rj k we find 
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where 



E k {iJ) 



7* 



7w7s 



(S14) 



(S15) 



Here £, C, and 7 stand for typical values. 

Actually, all impurities contribute to the deviation from inversion symmetry at site i, 
and therefore to 7 W . However, since this contribution depends on 1/Rf k , and has random 
sign, 7 W is strongly dominated by the contributions from the closest impurities. Thus we 
can assume that impurity k is a nearest neighbour impurity to impurity %. 
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(a) (b) 

FIG. S2: The true iS-spin density of states ng(E), as a function of energy for different short- 
distance behaviours of the interaction. Details are as in Figure 4 in the main text, only here x= 0.2 
and lattice size is 18 3 cells (ie., ~ 4650 TLSs). 

Let us now estimate the magnitude of 7 W . First, C ■ (5a/a) 2 ~ 5Muj 2 (Sa) 2 , where 5a 
denotes a typical strain, and the latter expression is the difference in kinetic energy of the 
impurity compared to the host ion. Approximating 5M ps M and uj ~ c/a we obtain 
C ~ Mc 2 . Since Rf 3 ~ a 4 n _4//3 , where n is the dimensionless impurity concentration, we 
get 7 W w n i/3 (/a. 

Let us now estimate (. Since ( is the coefficient of the second derivative of the dis- 
placement, it involves interaction with next nearest neighbors (see Fig. SI). Thus, 
C ■ 5a/a 2 ~ E$ ■ 5a/a so that £ w j s aE$/Ec- Thus, for n ~ 1 we find 7 W ~ •y s E^ > /Ec. 



1(3 



Note, that the ratio between the typical strain in a strongly disordered or amorphous solid 
to the interatomic distance, is 5a/a ~ E^/Ec, and therefore g = 7 w /7s ~ Sa/a. 

The calculation above can be similarly done for the effective random field exerted by the 
first term in the interaction Hamiltonian (IS4I) . with no reference to the S'-TLSs. 

1. dilute impurities in a lattice 

The reduction of x below x c ~ 0.2 has two main effects. First, the typical strength of 
the interaction of the r TLSs with the lattice becomes smaller, 7 W « g^ s n i ^. This, together 
with the dependence of the typical TLS-TLS interaction on n, would limit the universal 
regime to T < n 7 ^Tu- However, even within this much smaller temperature regime one 
can not expect universal behavior. This is since the variance of 7 W becomes large, and 
correlated with the TLS energy: strong local strains acting on the r spins from nearby 
S'-impurities increase both 7 W and J St . Such a correlation invalidates the assumptions of 
the ST model, and results in a temperature dependence of the phonon attenuation. In the 
mixed crystals (e.g. KBr:KCl with CN impurities) the situation is different, as the random 
fields experienced by the r TLSs are x independent. Thus, the condition required by the 
ST model that the impurities have similar interaction with the lattice, and that these are 
not correlated with their excitation energies, are satisfied. Furthermore, the disorder in 
the energies E T is dictated by the random fields, and is therefore independent of x. The 
only effect of having small x is the reduction in n(r), and we therefore find that universal 
properties exist in the mixed crystals, albeit with a tunneling strength C ~ 0.1gx/x c . 

C. Dipolar gap for the interacting S-r system 

In this section we calculate the gap for single particle excitations for both the S and r 
TLSs. Consider the Hamiltonian 



where = cf-J^ b /{R^ + a^), a,b = St, cfj 1 is a random number of order unity, J T Q T = 
gJg T = g 2 Jo S , and g <C 1 (g ~ 0.01 — 0.03 in our case). We introduce here a cutoff length 
a Q for the interaction. Rounding the diverging interaction at small distances is crucial for 




(S16) 
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the calculation of the gap of the S"-TLSs, but its functional form and the exact value of the 
cutoff parameter are only of small quantitative relevance. 

1. gap of t impurities 

Consider first the gap exerted by the S spins on the r spins, neglecting the r— r interaction 
(third term). Following Efros and Shklovskii, 



Here n°(E) ~ l/gJ corresponds to the ungapped DOS, and Uj is the interaction between 
the given r spin and the j'th S spin. We would like to show that there is no appreciable 
reduction of the r DOS, at least outside a window which is exponentially small in g. We 
proceed by assuming a series of simplifying conditions, all leading to a reduction of P(E T ) 
below its actual value, and obtain our final result in form of an inequality. 

We first overestimate ris(E) = 1/J Q for all E < J Q , neglecting their reduced DOS at low 
energies. We then enumerate the S spins according to their distance from the r impurity. 
Since there are j impurities within a volume (Vj) 3 , the maximum interaction is given by 
jjmax _ gj o jj_ We assume that all interactions have this maximum value, taking the short 
distance cutoff a a = and cfj = 1. Under the above assumptions, and defining E T = e, we 
obtain for e <C gJ Q 



n 



, T {E) = n° T (E) J] Q(E T + E Sj - 2U 3 ) = n°(E)P(E T ). 



(S17) 



j 
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Defining e' = e/2gJ Q we obtain 




(S19) 




Multiplying the denominators, and expanding in a series in g, we see that 
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(S20) 



and therefore 



P(e') > 1 -Ag\nl/e. 



(S21) 



i.e. P(e') « 1 for e' > exp (-1/8^). 

Below we will argue that the S impurities are strongly gapped themselves, thus also the 
above small correction at g 2 J Q < E T < gJ Q is probably an overestimate. 

To calculate the gap that the r impurities exert on themselves, one can repeat the same 
arguments as above, with s = J Q replaced by jf T = gJ Q and J^ T replaced by J T Q T = g 2 J Q . 
Eq. (IS2ip is thus reproduced as an inequality. However, since unlike the S impurities, the r 
impurities are not strongly gapped, one can recover the line of arguments with « instead of 
>, and conclude that 



where e" = e/g 2 J Q and c ~ g. Thus, the correction to the DOS at E < g 2 J a is reduced by 
the small parameter g, which is the ratio between the r — r interaction and their energy 
disorder, and is also the relevant small parameter in our theory. Experimentally, the dipole 
gap is indeed seen at energy scales comparable to the r — r interactions^, and its magnitude 
is considerably reduced. 

2. gap of S impurities 

The calculation of the gap of the S impurities is more subtle, since there is no small 
parameter, and the result is of order unity. Furthermore, the correlation between the in- 
teraction and r energies is crucial, and must be taken into account. We therefore solve the 
problem using a numerical simulation. Neglecting the third term in Eq. (1S16j) the values of 
E T are dictated by the second term. We are interested in P(E$), defined by 




(S22) 




(S23) 
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Here E is an energy a few times larger than Tjj and Uj is the interaction between a given 
S*-TLS and Tj. We therefore use the following algorithm: 

We choose a lattice size and concentration, and randomly distribute the impurities in the 
lattice. Each impurity has an S spin and a r spin. For a given a a , the interaction is given 
by Uij = CijS*Tj Jq t I '(Rfj + a„). We define = CijS*, where c^- is chosen from a random 
distribution with width unity. We then obtain and E Tj = £\ Uij with their essential 
dependence, but with no reference to the spin configuration of the S TLSs. We then flip the r 
spins where E T < 0, to have a positive excitation energy. The Uy's are accordingly redefined. 
For each Si we then calculate E^ x , the maximal Es { satisfying f| . Q(Es t + E Tj — 2Uij) > 0. 
We then obtain P(Es) = N(E^ X < Es)/N(S), the ratio between the number of spins with 
El? ax < Es to the total number of spins. The results for x= 0.5 are given in Fig. 4 in the 
main text, and for x= 0.2 in Fig. S2. 
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